Reionization of Hydrogen and Helium by Early Stars and Quasars 

J. Stuart B. Wyithe^'^ and Abraham Loeb^'^''* 
swyitheOcf a. harvard. edu; loebOsns . ias . edu 

ABSTRACT 

We compute the reionization histories of hydrogen and hehum due to the ionizing 
radiation fields produced by stars and quasars. For the quasars we use a model 
based on halo-merger rates that reproduces all known properties of the quasar 
luminosity function at high redshifts. The less constrained properties of the 
ionizing radiation produced by stars are modeled with two free parameters: (i) 
a transition redshift, z^^^^, above which the stellar population is dominated by 
massive, zero- met allicity stars and below which it is dominated by a Scalo mass 
function; (ii) the product of the escape fraction of stellar ionizing photons from 
their host galaxies and the star-formation efficiency, /esc/*- We constrain the 
allowed range of these free parameters at high redshifts based on the lack of the 
HI Gunn-Peterson trough at 2; < 6 and the upper limit on the total intergalactic 
optical depth for electron scattering, Tes < 0.18, from recent cosmic microwave 
background (CMB) experiments. We find that quasars ionize helium by a redshift 
z ~ 4, but cannot reionize hydrogen by themselves before ^ ~ 6. A major 
fraction of the allowed combinations of fescf* and Ztran lead to an early peak in 
the ionized fraction due to metal-free stars at high redshifts. This sometimes 
results in two reionization epochs, namely an early HII or Helll overlap phase 
followed by recombination and a second overlap phase. Even if early overlap is 
not achieved, the peak in the visibility function for scattering of the CMB often 
coincides with the early ionization phase rather than with the actual reionization 
epoch. Consequently, does not correspond directly to the reionization redshift. 
We generically find values of Tes > 7%, that should be detectable by the MAP 
satellite. 
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1. Introduction 

Following cosmological recombination at a redshift z ~ 10'^, the baryonic gas filling up 
the universe became predominantly neutral. Given that this gas is known to be mostly ion- 
ized today, one arrives at two of the major questions in current extragalactic astronomy: (%) 
when were the cosmic hydrogen and helium re-ionized? and (ii) which sources dominated this 
reionization process? The answers to both questions are likely to be different for hydrogen 
and helium (see review by Barkana & Loeb 2001). Recent observations provide preliminary 
answers to the first question. The absorption spectra of SDSS quasars at z ^ Q indicate that 
the neutral fraction of hydrogen increases significantly at z > 6 (Becker et al. 2001; Fan 
et al. 2002), and the UV spectrum of quasars implies that helium is fully ionized only at 
z < 3 (Jacobsen et al 1994; Tytler 1995; Davidsen et al. 1996; Hogan et al. 1997; Reimers 
et al. 1997; Heap et al. 2000; Kriss et al. 2001; Smette et al. 2002). The latter observations 
also indicate, through a cross-correlation between the hydrogen and helium forest of absorp- 
tion fines (Kriss et al. 2001; Smette et al. 2002), that the ionization of hefium at 2; ~ 3 
had significant contributions from both quasars and stars. No analogous evidence exists for 
hydrogen at 2; > 6. 

In this paper we make theoretical predictions for the reionization histories of hydrogen 
and helium. We calculate the contribution from quasars using a model that matches all the 
existing observational data on the quasar luminosity function at high redshifts (see Wyithe 
& Loeb 2002 for details). We model the less constrained stellar contribution using two 
free parameters, which are in turn constrained by existing observational data. Since the 
recombination times of both hydrogen and helium are shorter than the age of the universe 
in overdense regions of the intergalactic medium (IGM), it is possible that these species 
experienced more than one epoch of reionization. One of the goals of our detailed study is 
to identify the parameter values for which multiple reionization epochs are possible. 

The reionization history has important implications for the temperature and polar- 
ization anisotropics of the cosmic microwave background (CMB). Anisotropics on scales 
smaller than the angular size of the horizon at hydrogen reionization are suppressed by 
a factor ~ e""^"^ where Tes is the fine-of-sight optical depth for electron scattering (Hu & 
White 1997; Haiman & Loeb 1997), and secondary anisotropies are added (Hu 2000 and 
references therein). Reionization is the primary source of polarization anisotropies on large 
angular scales (Zaldarriaga &: Seljak 1997; Kamionkowski et al. 1997; Hu 2000). Polarization 
anisotropies might be detected in the near future by MAP on large angular scales (Kapling- 
hat 2002) or by ground-based experiments on small angular scales. So far, the available data 
on the temperature anisotropies of the CMB provides an upper limit of Tcs ^ 0.18 (Wang, 
Tegmark & Zaldarriaga 2002; Bond et al. 2002). We will use this upper limit to constrain 
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the free parameters of our model. 

Throughout the paper we assume density parameters values of flm = 0.35 in matter, 
Qft = 0.052 in baryons, Qa = 0.65 in a cosmological constant, and a Hubble constant of Hq — 
65 kms~^ Mpc~^ (or equivalently h — 0.65). For calculations of the Press-Schechter (1974) 
mass function (with the modification of Sheth & Tormen 1999) we assume a primordial 
power-spectrum with a power-law index n = 1 and the fitting formula to the exact transfer 
function of Cold Dark Matter, given by Bardeen et al. (1986). Unless otherwise noted we 
adopt an rms amplitude of erg = 0.87 for mass density fluctuations in a sphere of radius 



2. Reionization in a Clumpy Universe 

The simplest estimate of the epoch of reionization is based on the following simple 

considerations. Given a co-moving density of ionizing photons in a homogeneous but 
clumpy medium of comoving density no (where the size of the Hll region is much larger 
than the scale length of dumpiness), the evolution of the volume filling factor Qi of ionized 
regions is (Haiman & Locb 1997; Madau ct al. 1999; Barkana & Loeb 2001) 

dQi 1 duy C dt 

~T~ = 1 aB^Qifie-j-, 

dz Hq dz a"* dz 

where ckb is the case B recombination coefficient, a = 1/(1-1-2;) is the scale factor, rig is 
the comoving electron density, and C = {n^) / {uq)'^ is the clumping factor. This equation 
describes statistically the transition from a fully neutral universe to a fully ionized one, and 
yields reionization redshifts for hydrogen of around 7-12. Large uncertainties arise in both 
the source term and in the value of the clumping factor (because more rapid recombinations 
lead to a slower evolution of Qi). 

A more realistic description of reionization in a clumpy medium is provided by the model 
of Miralda-Escude et al. (2000). In what follows, we draw primarily from their prescription 
and refer the reader to the original paper for a detailed discussion of its motivations and as- 
sumptions. The model assumes that reionization progresses rapidly through islands of lower 
density prior to the overlap of individual cosmological ionized regions. Following overlap, 
the remaining regions of high density are gradually ionized. It is therefore hypothesized that 
at any time, regions with gas below some critical overdensity Ai = Pi/ {p) are ionized while 
regions of higher density are not. The assumption of homogeneity in equation (1) implies 
that the volume filling factor equals the mass filling factor. Therefore, within the model of 
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Miralda-Escude et al. (2000) we replace dQi/dz by dFy^i^ij/dz^ where 

FM(Ai)= r"dAPv(A)A (2) 
Jo 

is the fraction of mass in regions with overdensity below Ai, and Pv(A) is the volume 
weighted probability distribution for A. The product Q-JJ can be rewritten 

where -F(Ai) is the fraction of the volume with A < A^ The quantity Q;B(l+^)^We j^' (iAPv(A) A^ 
is therefore the number of recombinations per atom in the IGM per second. Note that the 
term analogous to the clumping factor is calculated from the volume weighted distribution. 
The mass fraction FM(Ai) (or equivalently Ai) therefore evolves according to the equation 

dFM(Ai) ^ 1 dn^ g^ ^^^i) ^ (4) 
dz Uq dz a? ^ dz 

This equation assumes that all ionizing photons are absorbed shortly after being emitted, so 
that there is no background ionizing field, and no loss of ionizing photons due to redshift. 
We therefore implicitly assume that the mean free path of ionizing photons is much smaller 
than the Hubble length. This should be valid at redshifts not too much smaller than the 
overlap redshift. 

The integration of equation (4) requires knowledge of Pv(A). Miralda-Escude et al. (2000) 
found that a good fit to the volume weighted probability distribution for the density as seen 
in N-body simulations has the functional form 

A-'^rfA, (5) 

with 5q = 7.61/(1 + z) and [3 = 2.23, 2.35 and 2.48, and Cq = 0.558, 0.599 and 0.611 at 
z — 2, 3 and 4. At 2; = 6 they assume /? = 2.5, which corresponds to the distribution of 
densities of an isothermal sphere, and solve for A and Cq by requiring the mass and volume 
to be normahzed to unity. We repeat this procedure to find Pv(A) at higher redshifts. The 
proportionality of Sq to the scale factor is expected for the growth of structure in an fl^ = 1 
universe or at high redshift otherwise, and its amplitude should depend on the amplitude 
of the power-spectrum. The simulations on which the distribution in Miralda-Escude et 
al. (2000) was based assumed Q^n = 0.4 in matter, Q\ = 0.6 in a cosmological constant and 
(Tg = 0.79, close to the values used in this paper. 



Pv(A)dA = Aexp 
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Equation (4) provides a good description of the evolution of the ionization fraction 
following the overlap of individual ionized bubbles, because the ionization fronts are exposed 
to the mean ionizing radiation field. However prior to overlap, the prescription is inadequate, 
due to the large fluctuations in the intensity of the ionizing radiation. A more accurate model 
to describe the evolution of the ionized volume prior to overlap was suggested by Miralda- 
Escude et al. (2000). In our notation the appropriate equation is 

dz = n^^- + ^(AcO-eQi^. (6) 

or 



dQi 1 dn. 



7 



dz nOFM(Acnt) dz 



aB(l + ^)^i?(Acrit)riei- + 



(7) 



-^M(Acrit) 



' dz dz 

In this expression, Qi is redefined to be the volume filling factor within which all matter 
at densities below Acrit has been ionized. Within this formalism, the epoch of overlap 
is precisely defined as the time when Qi reaches unity. However, we have only a single 
equation to describe the evolution of two independent quantities Qi and Fm- The relative 
growth of these depends on the luminosity function and spatial distribution of the sources. 
The appropriate value of Acrit is set by the mean free path of the ionizing photons. More 
numerous sources can attain overlap for smaller values of Acrit- Assuming Acrit to be constant 
with redshift, we find that results do not vary much (less than 10% in the optical depth to 
electron scattering) for values of Acrit ranging from a few to a few tens. At high redshift, 
these Acrit correspond to mean free paths comparable to the typical separations between 
galaxies or quasars. We assume Acrit = 20 (which lies between the values for galaxies and 
quasars) throughout the remainder of this paper. 



3. The Reionization of Hydrogen and Helium 

Next we describe the joint evolution of the filling factors and ionized mass fractions for 
hydrogen and helium using generalizations of equations (4) and (7). Let Qn+i Qne+ ^-nd 
QHe++ be the filling factors of ionized hydrogen, singly ionized helium and doubly-ionized 
helium, and Fh+ , -Fhc+ and -FHe++ be the mass-fractions of ionized hydrogen, singly ionized 
helium and doubly ionized helium, within the volumes Qh+j Qhc+ a-^d Qhc++ respectively. 
For later use we also define^ the fractions of ionized mass in the universe (5^+, Q^^+ and 



^Note that is the mass-fihing factor. The mass fiUing factors for H+ and He++ equal Qh+-Ph+ and 
<9He++-^He++ respectively. Prior to the overlap of Helll regions, Q^^+ = Q-ae+Pne+- However following the 
overlap of Helll regions, <9He+ ~ -^He+ • 
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We assign the comoving densities of hydrogen and hehum; — 1 .88 x 10'~^(Qb^^/ 0.022) cm~^ 
and = 0.19 x 10~^(f2b/?.^/0.022) cm""', respectively, and adopt the case B recombi- 
nation coefficients^ for H"*", He^ and He^+ at a temperature of lO^K (Osterbrock 1974; 
Barkana & Locb 2001); a^^ = 2.6 x lO'^^cm^ s-\ a^"^ = 2.73 x IQ-^-'^cm^ s^^ and 
Og^^^ = 13 X 10^^'^ cm^ s~^. The ionizing radiation field is described by the comoving den- 
sities of ionizing photons for ionized hydrogen, singly ionized helium and doubly ionized he- 
lium, namely n^^ , n^^^ and nj^°^^ , respectively. These quantities are computed through inte- 
gration of the corresponding spectra in the frequency intervals 3.29 x 10^^ < u < 5.94x lO^^Hz 
for HI, 5.94 X 10^^ <u< 1.31 x lO^^Hz for Hel (He+), and u > 1.31 x lO^^Hz for Hell (He++). 
Two classes of sources, quasars and stars, contribute to the ionizing radiation field. We dis- 
cuss these in turn in the following subsections. 



3.1. Ionizing Photons from Quasars 

To calculate the comoving density of ionizing photons emitted by quasars, n^, we inte- 
grate over both the quasar luminosity function and over the mean quasar spectrum (specified 
in Schirber & Bullock 2002). Given the B-band quasar luminosity function $(-Lb, z) at red- 
shift z in units of Mpc~^Lg^, we find for example 

— ^ = / dLB^LB,z) dv— , (8) 

az aZ Jo J3.29xl0i5Hz '^p'^ 

where is the Planck constant, and L^{Lb) is the luminosity in units of erg s~^ Hz^^ of 
a quasar with a B-band luminosity Lb- At redshifts below z = 2.5 we use the luminosity 
function of Boyle et al. (2000). At z > 2.5 we use the theoretical luminosity function 
of Wyithe & Loeb (2002); the reader is referred to the original paper for more details. 
This quasar evolution model successfully describes all known properties of the high redshift 
quasar luminosity function (Fan et al. 2001a,b), and reproduces measurements of the black- 
hole/dark matter halo mass relation (Ferrarese 2002) as well as estimates of the quasar duty 
cycle (Steidel et al. 2002). Since quasars at high redshift appear similar to typical examples 
at low redshift (Fan et al. 2001b), the luminosity function fully determines the ionizing 
radiation field. We therefore have almost no freedom in our calculation of the ionizing 
radiation field from quasars at high redshifts. This allows us to explore the effect of different 
stellar radiation fields on the reionization history of the universe. 



^At high redshift the absorbtion mean-free-path divided by the speed of hght is much smaller than the 
Hubble time. The appropriate recombination coefficient is therefore case B. 
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3.2. Ionizing Photons from Stcirs 

To compute the comoving densities of stellar ionizing photons, n^, we integrate over 
model spectra for stellar populations in the frequency ranges quoted above to find (for each 
species) the total number of ionizing photons, A^^, emitted per baryon incorporated into 
stars. We find 

= N^f^^J^—n^, (9) 

where rib is the local number density of baryons, /esc is the escape fraction for the ionizing 
radiation, f^, is the fraction of baryons within a galaxy that are incorporated in stars, and 
Fb(^) is the fraction of baryons in the universe at redshift z that collapsed and cooled inside 
galaxies. To compute F\)[z) we require the dark matter collapse fraction in halos Fco\{z, T-cain) 
with virial temperatures larger than Tmin, given by 

1 f°° dn 
F.oi(z, Tlnin) = — / dMM^, (10) 

Pm J M ■ (T ■ ) dm 

where pm is the co-moving mass density in the universe, and drips/ dM is the Press-Schechter (1974) 
mass function for dark matter halos. In regions of the universe that are neutral, the 
critical virial temperature (the temperature at which atomic cooling becomes efficient) is 
Train — lO^K (a circular velocity of fcir ~ 10km s"-*^). In reionized regions, infall of gas is 
suppressed below a virial temperature of Tmin ~ 2.5 x lO^K [vdr ~ 50km s~^, Thoul & Wein- 
berg 1996). Given Fco\{z, lO^K) and Fco\{z,2.b x l^fK) we can take the baryonic collapse 
fraction 

Fb(z) = QS+Feoi(-2, 2.5 X lO^K) + (1 - Q^+)F,,i{z, lO^K). (11) 

This prescription for calculating the stellar ionizing field will be referred to hereafter as case 
A. Note that /esc could be different for hydrogen and helium; however, we adopt a single 
value for it assuming that this value is dictated by the geometry of the optically-thick gas 
that allows limited escape routes which are common for all ionizing photons (Wood & Loeb 
2000 and references therein). 

We also consider a second evolution for the stellar ionizing field, denoted hereafter as 
case B. A recent study (Kauffmann et al. 2002) has shown that in a large sample of local 
galaxies, the ratio e = M^/Afhaio (where and Mhaio are the total stellar and dark matter 
halo masses respectively) scales as e oc M^^J, for < 3 X IO^^Mq, but is constant for 
larger stellar masses. Note that the star-formation efficiency is proportional to e. Since 
star-formation is thought to be regulated by supernova feedback (Dekel & Silk 1986), the 
important quantity is the depth of the galactic potential well, or equivalently the halo circular 
velocity. Using the stellar mass TuUy-Fisher relation of Bell & De Jong (2001), we find the 
threshold circular velocity v^, — 176 km s~^ that at z = corresponds to a stellar mass 
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of 3 X 1O^°M0. In this case we define fi, as the star- formation efficiency in galaxies with 
circular velocities larger than v^, and calculate the effective product of the dark matter 
collapse fraction and star-formation efficiency 



1 r°° 



where e 
calculated from 



1 for M > M^i^^ and e 



(M/M4iJ^/^ for M < MU- The value of M^^^J 



(12) 



z) is 



4.3 X 10^°/i 



V 176km s" 



crit 



-1/2 



l + Z 

10 



-3/2 



Mo, (13) 



where A^it = 1877^ + 82d - 39c^^ d = - 1 and = + z)V[Qin(l + -z)^ + ^^a]- 
Finally, the value of /*Fb required for calculation of dn^/dz may then be computed in 
analogy with case A, and substituted into equation (9). 

In both cases A and B we make the distinction between the ionizing radiation field 
due to a possible early population of zero-metallicity stars, and the metal enriched stars 
observed at lower redshifts. It is thought that the primordial initial mass function favored 
massive stars (Bromm, Copi & Larson 1999, 2001; Abel, Bryan & Norman 2000; Mackcy, 
Bromm, Hcrnquist 2002). The possible existence of this population is very important for 
reionization because the spectrum of these stars would result in an order of magnitude more 
ionizing photons per baryon incorporated into stars (Bromm, Kudritzki, & Loeb 2001). The 
formation of the very massive stars was suppressed as the material out of which stars form 
was enriched with metals. The fraction of the ionizing photons produced by metal-free 
stars depends on several unknown parameters, including the mixing efficiency of metals, 
the environments in which new stars form, and most importantly, the threshold metallicity 
above which star formation is transformed from being dominated by massive stars to a 
Scalo (1998) initial mass function (IMF). The threshold metallicity is believed to be small; 
Bromm et al. (2001) argue for a threshold < 10^^ of the solar metallicity value. The 

efficiency of mixing of metal enriched outflows from star forming galaxies to the surrounding 
IGM is even more uncertain; Scannapieco, Ferrara & Madau (2002) find that the mass 
weighted mean metaUicity can reach values greater than 10~^ of the solar value at redshifts 
as high as 20, and note that the average metaUicity scales with star formation efficiency, 
supernovae rate, and the fraction of supernovae energy that is channeled into outflows. The 
average metallicity increases roughly exponentially with redshift as it is modulated by the 
exponential growth in the collapse fraction of baryons at high redshifts. 



Since the first metals were produced by supernovae, it is reasonable to suppose that 
the enrichment of the IGM with metals at redshift z is proportional to the mass of baryons 
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that has formed stars by that redshift. The average metaUicity at a redshift z can then be 
written as 

^ = C^eta.%^ Tci.^, (14) 

Zjq Zjq az 

where Cmetai is a constant. Note that the baryonic coUapse fraction Fb, which is computed 
from equation (11) is affected by the reionization history. The formation of metal-rich stars 
is assumed to be proportional to the average metalicity, with the fraction reaching unity 
when the average metalicity of the IGM reaches Before reaches ^^^^ we may 

approximate the fraction of star formation at redshift z that is in zero-metallicity stars to 
be (1 — :^£l / ^thresh ') ■vv'hile after reaches ^i^^ we may assume all star formation to have 
a Scalo (1998) IMF^. 

For a given star formation efficiency and escape fraction, the above prescription yields a 
characteristic redshift for the build-up of the threshold metal enrichment. The exponential 
growth in the collapse fraction with redshift implies that the evolution in the mode of star 
formation resembles a step- function. We therefore define a transition redshift ^tran below 
which metal-rich stars with a Scalo (1998) IMF dominate the production rate of ionizing 
photons. 

Through most of the paper we present results in terms of -Ztj-anj but return to a justifi- 
cation of this choice in § 7. At redshifts above -Ztraru we use values of A^^ calculated for each 
of the three species H"*", He"*" and He"'"^ assuming massive (> lOOM©) zero-metallicity stars 
and the generic spectrum calculated by Bromm, Kudritzki & Loeb (2001). The resulting 
values of are ~ 14020 for H+, ~ 25200 for He+ and ~ 4480 for He++. At lower redshifts 
we assume metal enriched stars (l/20th solar metaUicity) with a Scalo (1998) mass-function, 
and use spectral information from the stellar population model of Leitherer et al. (1999)^. 
This results in values of A^^ ~ 3250 for H+, A^^ ~ 1020 for He+ and A^^ ~ 0.2 for He++. 



3.3. Constraints on the Evolution of H"*", He"*" and He"'"+ Ionization Fronts 

The typical distance that an ionizing photon of neutral helium can travel through neutral 
hydrogen having a cross-section a-^{y) before the optical depth reaches unity is the absorption 

''The situation in reality is more complicated since mixing of metals is incomplete and the formation 
sites of new stars are correlated with the enriched regions. Our discussion uses the minimum number of free 
parameters to describe this complicated process. 

^Model spectra of star-forming galaxies were obtained from http://www.stsci.edu/science/starburst99/. 



-10- 



Case I 



Case II 



high density {A>i 


^crit) 








Acrit 










m 

source 










• 

source 










^crit 






high density (A> 


Acrit) 





source 



^crit 




^crit 





source 




□ dense regions; 
neutral H, neutral He 

□ neutral He 

H+ He+ 
^ H+He++ 
# Ionizing sources 



Fig. 1. — Schematic diagram of the process of reionization. Prior to the overlap of hydrogen, 
reionization proceeds through low density regions (below Acrit) for both hydrogen and helium. 
Following the overlap of HII regions, the hydrogen reionization front moves into denser 
regions but helium reionization continues to proceed through regions with density below 
Acrit- Finally, after the overlap of Helll regions, the reionization of helium is free to proceed 
into denser parts of the IGM. 



mean-free-path A, where 





Wc find A ~ 4 [(1 + z) /lO] 



the typical size of an HII region around a IO^Mq halo prior to reionization (Barkana & 
Loeb 2000). Therefore, at z ~ 10 the He"*" propagation front can only lead the H"*" propa- 
gation front by a small fraction of the size of the HII region, and the excess He"*" ionizing 
photons photo-ionize additional hydrogen. On the other hand, most ionizing photons for 
hydrogen (z/ > 5.94 x lO^^Hz) are not absorbed by hehum, and so the expansion of the HII 
region is not inhibited by the presence of helium. Obviously, due to the absence of singly 
ionized helium the He++ front cannot propagate beyond the H+ front. 

Following these considerations, we impose two restrictions on the evolution of (5He+ 
QHe++- First, Q}ie++ ^ Qh+- Second, QHe+ < Qh+ — QHe++- Following helium overlap, we 
impose the same restrictions on the evolution of the mass fraction, i.e. FHe++ < i^H+ and 
-^He+ < -^H+ — -^He++ ■ The fuU equations governing the evolution of the filling factors and mass 
fractions are presented in the Appendix. The equations are split into three cases describing 
three different epochs. The evolution of the filling factors prior to hydrogen overlap are 
described in case I. Case II describes the evolution of the filling factors and of the ionized 
hydrogen mass fraction following the overlap of HII regions, but prior to the overlap of Hell 
regions. Finally, case HI describes the evolution of the ionized mass-fractions following the 
overlap of cosmological Hell regions. Figure 1 provides a schematic representation of the 
three epochs in reionization history. 



In this section we derive reionization histories for different values of Znan and /esc/*- 
These histories are shown in figures 2 (for case A star-formation) and 3 (for case B star- 
formation) . 

For case A wc show results for combinations of /esc/* = 0.0005 and 0.005, and 2;t,ian = 
10.0 and 20.0. Wc note the following features. Large values of /esc/* allow the massive zero 
metallicity stars to reionize both hydrogen and helium in the universe at very early times. If 
the zero metallicity stars do not reionize helium, quasars result in the overlap of Hell regions 
by 2; ~ 4, consistent with observations showing transmission just blueward of the helium 
Lya line at 2; ~ 3 (Jacobsen et al 1994; Tytler 1995; Davidsen et al. 1996; Hogan et al. 
1997; Reimers et al. 1997; Heap et al. 2000; Kriss et al. 2001; Smette et al. 2002) as well as 
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Reionization Histories 
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Fig. 2. — Sample reionization histories assuming case A star formation. Dark lines represent 
the evolution of the volume filling factors (5h+ , QHe+ and Q}ie++ , while the light lines show 
the evolution of the fractions of ionized mass in the universe Q^+, Q'^g+ and Q^^++- The 
solid, dashed and dotted hues refer to H+, He+ and He++, respectively. 
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Fig. 3. — Sample reionization histories assuming case B star formation. Dark lines represent 
the evolution of the volume filling factors (5h+ , QHe+ and Q}ie++ > while the light lines show 
the evolution of the fractions of ionized mass in the universe Q^+, Q'^g+ and Q^^++- The 
solid, dashed and dotted hues refer to H+, He+ and He++, respectively. 
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evidence for a temperature rise in the IGM at z ~ 3.4 indicating the reionization of hehum 
at that time (Schaye et al. 2000; Theuns et al. 2002). Interestingly, we also find that it is 
possible for hydrogen and/or helium to have been reionized twice. This situation arises if the 
zero-metallicity stars reionize the universe prior to -Ztj-am but below Ztran the recombination 
rate is sufficiently high that the softer spectra of normal stars cannot maintain overlap. 
Overlap of HII regions is then re-achieved at a later time as the collapse fraction grows and 
the density of the IGM is lowered. A second overlap of Hell regions results from the rise in 
quasar activity at z < 10. Similar results are found in case B, for values of /esc/* that are a 
factor of ~ 10 larger than in case A. 

We would like to reiterate the point that multiple reionization epochs occur for some 
sets of parameters as a result of multiple peaks in the emissivity of ionizing photons as a 
function of cosmic time. In such cases the first emissivity peak occurs at -Ztran, at which 
point it drops due to the softer spectra of metal enriched stars. Later, when the collapsed 
fraction of baryons has reached a sufficiently high value, the resulting increase in emissivity 
further reionizes the IGM. As we will show in § 7, this double peaked behavior is still 
seen if a smoothly varying fraction of collapsed objects at mctallicities below the threshold 
produce the hard spectra rather than sources emitting prior to a fixed -Ztran- A double peaked 
emissivity will always arise so long as the relevant metal enrichment time is less than the 
time required for the collapse fraction to change by a factor equal to the ratio between the 
emissivities of the first generation of stars and subsequent metal enriched stars. ^ 

There arc two small, sharp jumps in the reionization histories that arc artifacts, but 
warrant explanation. The first is in the evolution of Q^^^ immediately following the overlap 
of Helll regions. This small jump (note that the curve is smooth, not discontinuous) arises 
because our formalism demands (5^6+ = at this epoch. However following the overlap of 
He"^"*", the He"*" front is free to rapidly expand, powered by the strong ionizing radiation field, 
and the mass fraction FHe+ therefore grows until the He+ front reaches the H"*" front. The 
second jump occurs at ^ = 2.5 where we switch from computation of the ionizing fiux due 
to quasars using the model quasar luminosity function, to calculation using the empirical 
luminosity function. Neither is a perfect description around this redshift, and there is a 
discontinuous jump in the luminosity density, and hence in the ionizing radiation density, 
resulting in a jump in the reionization histories. 

At early times the rise in temperature associated with reionization curtails star-formation. 



''Note that even if the mixing of metals is partial and the intergalactic filling factor of metals is smaller 
than unity (e.g., Scannapieco et al. 2002), new stars are likely to form in regions that were already enriched 
because galaxies tend to be clustered on large scale sheets and filaments. This effect needs to be included 
when evaluating the effective ^;tran from numerical simulations. 



- 15 - 



The effect of this is clearly seen in panel (a) of figure 2. Here the zero metaUicity stars rapidly 
increase the ionized fraction, but as Qh+ approaches unity, the ionizing flux is reduced, and 
the rate of reionization slows. This feature is not seen in cases where overlap occurs at lower 
redshift (e.g. panel (d) in figure 2). Many combinations of the parameters 2;tran and /esc/* 
predict significant reionization without overlap due to the early presence of metal-free stars, 
resulting in a peak in the ionized fraction centered on the transition redshift, followed by 
recombinations before further reionization due to normal stars and quasars. We term these 
peaks failed overlaps. In section § 6 we show that these failed overlaps leave a significant 
imprint on the visibility function for electron scattering of CMB photons. 



4.1. The Optical Depth to Electron Scattering 

Different reionization histories result in different densities of electrons as a function 
of redshift. One simple probe of the reionization history is the optical depth to electron 
scattering Tes, which depend on the mass fiUing factors Q'^ 

/■^°°° cdt 

Tes = dz—ar [Q^+nl + g^?,+< + 2QS++ny (1 + z)\ (16) 

where 0"^ = 6.652 x 10^^^ cm^ is the Thomson cross-section. In the following section we plot 
contours of Tes as a function of 2;tran and /esc/*- 



5. Joint Constraints on ztran and /esc/* 

In this section we discuss the constraints that current observations place on Zt^an and 
/esc/*- First we discuss the variation of Tes with Ztran and /esc/*- The optical depth for electron 
scattering is constrained from recent CMB anisotropy experiments to a value Tes < 0.18 
(Wang & Tegmark 2002). An upper limit on t^s implies that there is maximum redshift 
beyond which the universe was neutral. As a second constraint, we note that the spectra of 
quasars at z < 6 do not show a Gunn-Peterson (1965) trough, indicating that HII regions 
had achieved overlap earlier than that time. 

As we have seen in the previous section, the universe is reionizcd earlier for higher values 
of /esc/*; which rcsult in higher values of Tes. Hence /esc/* cannot be too large so as not to 
violate the fimit Tes < 0.18. On the other hand, if the product /esc/* is too small, then 
the ionizing radiation field will be too weak for the HII regions to overlap by 2; ~ 6. If the 
transition redshift 2;tran is large, then /esc/* may obtain higher values since the harder spectra 
of metal free stars are only available to the ionizing radiation field at high redshifts when 
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Fig. 4. — Contours of optical depth for Thomson scattering, Tes (sohd hnes) and contours 
of (5h+ at 2; = 6 (dashed hnes), as a function of /esc/* and -Ztran- Case A star- formation 
was assumed, and the shaded regions are excluded by existing observations. The 4 labeled 
points refer to the locations of the parameter sets used in the example histories presented in 
figure 2. 
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Fig. 5. — Contours of optical depth for Thomson scattering, Tes (sohd hnes) and contours 
of Qii+ at 2; = 6 (dashed hnes), as a function of /esc/* and -Ztran- Case B star-formation 
was assumed, and the shaded regions are excluded by existing observations. The 4 labeled 
points refer to the locations of the parameter sets used in the example histories presented in 
figure 3. 
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the IGM is dense and the collapse fraction is small. However if ^tran is low, then /esc/* may 
take smaller values and still result in overlap of HII regions by z ~ 6. In addition to the 
above constraints, observations suggest that the IGM was metal enriched by 2; ~ 6 (Songaila 
2001) to a level that would lead to a Scalo stellar mass function (Bromm et al. 2001) if new 
galaxies formed out of it. As a result, ^tran is limited to be larger than six. 

Figures 4 and 5 show contours of t^s (solid lines) as a function of Ztran and /esc/* for 
case A and B star-formation respectively. Also shown are contours of Qh+ at 2; = 6 (dashed 
lines), and a vertical line at ^tran = 6 (dot-dashed line). The regions excluded by the 
aforementioned constraints are shaded grey. For our case A star-formation model wc find 
that /esc/* < 0.004 if Ztri^n < 10. For a higher Zuan, we find that larger values of /esc/* 
are possible. If Ztrg,n ~ 6, /esc/* may be as low as 10"^, though for Zi^^^n > 15 we find 
/esc/* ^ 0.001. We find that the lower fimits on /esc/* are larger by a factor of ~ 3 for 
our case B star-formation model. However since the intensity of the stellar radiation field 
is curtailed by feedback at high redshift in this model, we find that /esc/* can be as large 
as ~ 0.1 for Ztran ^ 15, or even larger for Ztran ^ 15. Interestingly, the contours for Tes and 
^reion are uot parallel. A significant range of Tgs exist for different Ztra.n (particularly where 
case A star-formation is assumed). This indicates that measurement of optical depth does 
not determine the overlap redshift. We discuss the reason for this puzzling result in § 6. 

In § 4 we mentioned that hydrogen and/or helium could be rcionizcd twice, and showed 
examples of rcionization histories that exhibit this behavior. In figures 6 and 7 wc show 
contours of the overlap redshift for hydrogen (left panels), and helium (right panels), for case 
A and case B star-formation models. The redshift of the most recent overlap is described 
by the dark contours. At high values of Ztran and low values of /esc/*, the overlap redshift 
is almost independent of star-formation since then quasars dominate the ionizing radiation 
field. This is particularly true for helium, which is ionized by quasars at 2; ~ 4 in the 
absence of any stellar ionizing radiation. Cases where there was an additional earlier overlap 
are denoted by the grey contours. The region of the [^tran, /esc/*] plane that results in multiple 
overlap epochs is shaded grey. We see that this area covers a significant portion of parameter 
space, and therefore that it is quite possible the universe was reionized twice. The contours 
representing the upper limit of Tes < 0.18, and the lower limit for the overlap redshift z > 6 
are also shown (thick dashed lines). As mentioned earfier, parameters in the region between 
these limits are not excluded by current observations. 
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Fig. 6. — Contours of overlap rcdshift for hydrogen (left panel) and fully-ionized helium 
(right panel), as a function of /esc/* and ^tran assuming case A star-formation. The dark 
contours show the most recent overlap. Parameter pairs that produce two overlap epochs 
lie in the shaded region. The grey contours describe the redshifts of these earlier overlaps. 
Note that the regions below and to the right of the lowest redshift contour shown have nearly 
constant overlap redshift. This arises because quasars dominate the reionization process for 
these parameter pairs. Also shown for reference are the contours of Tes = 0.18 and (5h+ = 1 
at z — 6 (thick dashed lines) . Parameter pairs not lying between these lines are excluded by 
current observation. 
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Fig. 7. — Contours of overlap redshift for hydrogen (left panel) and fully-ionized helium 
(right panel), as a function of /esc/* and 2;tian assuming case B star-formation. The dark 
contours show the most recent overlap. Parameter pairs that produce two overlap epochs 
lie in the shaded region. The grey contours describe the redshifts of these earlier overlaps. 
Note that the regions below and to the right of the lowest redshift contour shown have nearly 
constant overlap redshift. This arises because quasars dominate the reionization process for 
these parameter pairs. Also shown for reference are the contours of Tes = 0.18 and Qii+ — 1 
at z — 6 (thick dashed lines) . Parameter pairs not lying between these lines are excluded by 
current observation. 
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6. The Visibility Function 

The visibility function g{z)dri describes the probabihty that an observed photon was 
scattered between conformal times r]{z) and rj + Arj 

dr 

9{z) = -^e-^- = aT [g^?+n° + +n°, + 2Q- ++n°,] (1 + zfe-^-, (17) 

where drj — {1 + z)cdt, Tes is the optical depth from redshift to z and the are mass 
filling factors. It is commonly assumed that the visibility function is peaked during the 
recombination era, with a second, lower amplitude, but broader peak coinciding with the 
epoch of overlap of HIT regions. However we have found that overlap may have occurred on 
multiple occasions, and that many combinations of the parameters ^rtran and fcscf* predict 
significant reionization without overlap {failed overlap) at high redshift due to the early 
presence of metal-free stars. During these failed overlaps, examples of which are shown in 
figures 2 and 3, the ionization fraction can be significant. Following a failed overlap at 
Ztran, recombinatious dominate the evolution and the ionization fraction drops before further 
reionization ensues due to normal stars and quasars. The higher cosmic density at early 
times might more than compensate for the smaller ionization fraction, so that the peak of 
the visibility function does not coincide with the epoch of reionization. 

This is indeed what we find. Figures 7 and 8 show visibility functions corresponding to 
the sample reionization histories shown in figures 2 and 3. Generally, if the early metal- free 
stars succeed in achieving reionization, then the peak in the visibility function coincides with 
this overlap epoch (e.g. case A in figures 2 and 3). However the optical depth Tcs is large and 
many of these cases arc already excluded by CMB measurements. Of the 8 cases shown in 
figures 2 and 3, those that exhibit failed overlap have visibility peaks that coincide with this 
failed overlap, rather than with the actual reionization epoch. Thus, the visibility function 
may probe the nature of the early generation of stars rather than the reionization epoch itself 

The reionizcd gas during a failed overlap will be in bubbles. Gruzinov & Hu (1998) have 
pointed out that a universe which is ionized in patches induces additional CMB anisotropics 
on arc-minute scales. Thus, patchy reionization due to a failed overlap at high redshifts will 
modify the shape of the small-scale power spectrum of GMB anisotropics. However, the 
amplitude of the added secondary anisotropies was predicted to be at a level well below the 
feature in the spectrum detected more recently by the Gosmic Background Imager (Padin et 
al 2001). This feature was explained by Bond et al. (2002) as being due to the foreground 
Sunyaev-Zeldovich effect from X-ray clusters. 
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Fig. 8. — Sample visibility functions (in units of Hq/c). The four cases correspond to the 
sample reionization histories (assuming case A star-formation) plotted in figure 2. 
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Fig. 9. — Sample visibility functions (in units of Hq/c). The four cases correspond to the 
sample reionization histories (assuming case B star-formation) plotted in figure 3. 
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7. Justification for Step-Function Approach to the Evolution in the Mode of 

Stcir Formation 

As mentioned in § 3.2 we have approximated the evolution in the mode of star formation 
as a step function. That is, we have assumed that the production rate of ionizing photons 
in the universe is transformed sharply at 2;tran from being dominated by massive, metal-free 
stars to being dominated by metal-rich stars with a Scalo IMF. In this section we justify this 
approach by comparing its results to those obtained with a scheme involving gradual metal 
enrichment of the IGM as described in § 3.2 by equation (14). 

Figure 10 shows two sample reionization histories computed using the gradual enrich- 
ment scheme, together with the resulting visibility function. The examples correspond to 
panels a) and c) for Case B star formation in Figure 3. The thick grey hne in these panels 
shows the fraction of star formation resulting in a Scalo (1998) mass function as a function 
of redshift. In each of these cases the value of Cmctai (listed in the figure) produces around 
50% Scalo (1998) IMF star formation at the previously used redshift -Ztran- Note the different 
values of Cmetab ciud the different histories of the collapsed mass- fraction. The features dis- 
cussed for the corresponding examples in Figures 3 and 9 are still present in Figure 10. These 
include the double overlap of He III in a) , and the non-coincidence of the overlap redshift of 
HII and the peak in the visibility function in b). Furthermore, we find unchanged values of 
Teg. These calculations appear to justify our use of a single parameter Ztran to characterize 
the reionization histories and visibility functions. 



8. Conclusion 

We have explored the reionization histories of hydrogen and helium due to stars and 
quasars. The results were analyzed as a function of the two free parameters in our model, 
namely: (i) the transition redshift, -^tran, above which the stellar population is dominated by 
massive, zero metallicity stars; and (ii) the product of the escape fraction of ionizing photons 
and the star-formation efficiency, /esc/*- The quasar model was not varied since it provides 
an excellent fit to all existing data on the luminosity function of quasars up to redshift ^ ~ 6 
(Wyithe & Loeb 2002). 

Figures 2 and 3 show sample reionization histories for different choices of Zt^an and 
/esc/*- We find that a wide range of 2;tran ^ 6 is allowed for /cscA ~ 2 x 10^^ if the 
star-formation efficiency does not depend on galaxy mass (figure 4) or /esc/* ~ 10^^ if 
the efficiency is suppressed in low mass galaxies (figure 5). This wide range satisfies the 
constraints that overlap of the HII regions must be achieved by ^ ~ 6 (Fan et al 2002) and 
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Fig. 10. — Results for reionization histories computed assuming continuous enrichment of 
the IGM rather than sudden enrichment at ^tran- The results should be compared with 
panels a) and c) of Figures 3 and 9. Left: Sample reionization histories assuming case B 
star formation. Dark lines represent the evolution of the volume filling factors Qn+, Qnc+ 
and Quc++, while the light lines show the evolution of the fractions of ionized mass in the 
universe (5h+' ^hc+ ^He++- '^^^ solid, dashed and dotted lines refer to H+, He"*" and 
He''""'", respectively. Right: The corresponding visibility functions (in units of Hq/c). 
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that the optical depth for electron scattering must be limited to Tes < 0.18. The allowed 
range leads generically to Tes ^ 7% (see figures 4 and 5). The MAP satellite is expected to 
have sufficient sensitivity to detect values of Tes ^ small as ~ 5% (Kaplinghat et al. 2002), 
well below this range of expected values. 

A major fraction of the allowed range of /esc/* and Ztran leads to an early peak in the 
ionized fraction due to metal-free stars at high redshifts. Often this peak results in a small 
but significant filling factor that is subsequently reduced (temporarily) due to recombination. 
In a restricted range of the allowed parameter values, we find that cither hydrogen or helium 
experience two overlap epochs, separated by recombination (sec shaded regions in figures 6 
and 8). That helium might have been reionized twice due to the presence of early very 
massive stars was previously pointed out by Oh et al. (2001). The first overlap phase is 
caused by the population of zero-metallicity, massive stars and the second is dominated by 
the quasars for helium or by stars and quasars for hydrogen. Even if early overlap is not 
achieved, the peak in the visibility function for scattering of the CMB often coincides with 
the early partial ionization peak rather than with the actual reionization epoch (see figures 7 
and 9). The resulting value of Tes is therefore larger than expected based the reionization 
redshift alone. Thus, the CMB visibility function may be probing the nature of the early 
generation of stars rather than the reionization epoch itself. 

Future CMB experiments, such as MAP and Planck, will provide tighter constraints 
on and will reduce the range of allowed /esc/* in om model (see figures 4 and 5). Any 
additional observational data on the composition of the stellar population or the abundance 
of Hel or Hell at z > 4 will provide stronger lower limits on Zuan- 
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Equations Describing Reionization of Hydrogen and Helium in a Clumpy 

Universe 



The filling factors Q}i+, Qhc+ and Quc++j and ionized mass fractions Fh+, -Fhc+ and 
-^He++ evolve according to six coupled ODEs. It is convenient to write the ODEs and con- 
straints in terms of the following quantities: 
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(Al) 



In the above definitions, the Sc and Scf are source terms describing the ionizing radiation 
field, while the Sk and Skf arc sink terms describing the recombination rate. In the recom- 
bination terms, the ric are the electron densities available for recombination (equal to the 
mean cosmological density multiplied by the volume filling factor and the mass fraction). 
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These are therefore 



H+ 



n. 



n: 



e 

and 

He++ 



, QHe+-^He+ + 2(5He++-^He++ 



(A2) 



We write the generahzation of equations (4) and (7) for the three species in terms of the 
above quantities. There are 3 cases corresponding to the period before HII regions overlap, 
the period following the overlap of HII regions but preceding the overlap of Hell regions, 
and the period following the overlap of Hell regions: 

•CASE I) Pre Hydrogen Overlap: If Qh+ < 1 or Fh+ < FM(Acrit) 

Prior to the overlap of HII regions, the mass fractions F and their derivatives are set 
by the critical value of Acrit chosen as the ionization density threshold. The fiUing factors 
Qn+ , Q}ie+ and Q}ie++ evolve according to a generalization of equation (4) . Hence 
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In the above, A 



dn: 



He+ 



and A 



are the excess ionizing photon rates having frequencies 



dz dz 

5.94 X 10^^ <v < 1.31 X lO^^Hz for Hel ionization, and v > 1.31 x lO^^Hz for Hell ionization. 
The values of these excess source terms are determined from the constraints on the relative 
evolution of the different ionization fronts. As these are excess rates, their values must be 
positive. 

If <5He++ < Qh+ then (5He+ is limited to be smaller than (5h+ — Qne++ since the He"*" 
front cannot propagate beyond the H+ front. Therefore, if (5He+ — Qn+ ~ QHe++ we require 
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in this case we find 



In addition to the above constraint, the He" 
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filhng factor Quc++ is hmited to be less than 
Qu+- Therefore, if Q}ie++ = Qu+ we require '^'^^"^^ = so that the fronts propagate 

at the same rate, and '^'^j^""'" — since the He"*" front cannot propagate if the Hell region is 
filhng all of the available volume. These constraints yield 
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•CASE II) Post Hydrogen Overlap: If Qh+ 

or FHe++ < -FM(Acrit)). 



1 and Fh+ > FM(Acrit) and {QHe++ < 1 



Following the overlap of HH regions, the mass fraction of hydrogen Fh+ is free to evolve 
according to equation (7) (with modified source term to account for the excess ionizing 
photons). The filling factors Qh+i Qhc+ and (5hc++ again evolve according to a generalization 
of equation (4). Note that if Qi = 1, substitution of equation (4) into equation (7) yields 
^ = 0. Thus, post HH overlap the value of (5h+ remains at unity. However by continuing 
to follow all the equations we allow for the possibility of a recombination epoch following an 
early reionization. We have 
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As before, the values of the excess ionizing photon rates are determined from constraints 
on the co-evolution of the ionization fronts. In the case where QHe++ < Qh+, Qhc+ is 
limited to be smaller than (5He+ = Qk+ ~ Qnc++- However as mentioned above, if (5h+ = 1 
then the above equations imply '^9m± — q. We therefore have ^^^^ — — ^5h£±+ ^ and since 



dz 



dz 



in this case we find 
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(AT) 



•CASE III) Post Helium Overlap: If Qhc++ = 1 and Fhc++ > Fm(A 



crit ) 



Finally, following the overlap of Hell regions, the mass fractions Fh+, i^Hc+ and -Fhc++ 
arc free to evolve according to equation (7) (with modified source terms to account for the 
excess ionizing photons). The filling factors Qh+i Qhc+ and Qnc++ still evolve according to 
generalizations of equation (4). At the point of He^+ overlap our formalism demands that 
there is no He"^ in the universe since Qne+ — 0. However, following He"^"*" overlap the He"^ 
front (which must lead the He^''" front) is free to propagate into the denser IGM. Therefore, 
following the overlap of He+"^, FHe+ becomes the fraction of helium in the universe (equal to 
zero initially) that is singly ionized. The recombination rate B^^^ is calculated 



AHe++ dA 



(A8) 



where Fhc++ = /(f ' d^^^ and Fhc++ + -FHe+ = /q 
substitution of equation (4) into equation (7) yields ^ 
(5He++ therefore remain at unity. We have 



dA^A. As before, if Qi = 1 
0. The filling factors (5h+ and 
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The excess ionizing photon rates are determined as follows. If FHe++ < then FHe+ is 
limited to be smaller than -Fhc+ = i^H+ — -^Hc++ since the He"*" front cannot propagate beyond 
the H"*" front. Therefore, if FHe+ > -Fh+ — -^He++ we require ^^^^ = '^^^^ — ^^^7-*-, and since 

= in this case, we find 
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f-^ + 4) ' \{S^f - SS) - {Sif - S^;)] . (AlO) 



The He mass fraction -Fhc++ is limited to be less than F^+. Therefore, if Fhc^ 
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so that the fronts propagate at the same rate, and 
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= Fh+ we 
since the 



He"*" front cannot grow if the He"*""*" is filling all of the available volume. These constraints 
yield 
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(All) 



